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This  paper  investigates  the  relation  between  microstructure,  macroscopic  transport  properties,  and  fab¬ 
rication  processing  for  a  gradient  porous  cathode  of  solid  oxide  fuel  cells  (SOFCs).  Functionally  graded 
porous  cathode  with  smooth  variations  in  pore  size  is  composed  of  lanthanum  strontium  manganite 
(LSM)  fabricated  on  yttria  stabilized  zirconia  (YSZ)  electrolyte  substrate  using  a  multi-step  spray  pyroly¬ 
sis  (SP)  technique  at  various  deposition  conditions.  Two-dimensional  (2D)  serial-sections  of  the  gradient 
porous  microstructure  obtained  by  FIB-SEM  are  fully  characterized  using  statistical  correlation  functions. 
Results  of  statistical  analysis  of  the  microstructures  revealed  that  the  SP  processing  technique  is  capa¬ 
ble  of  generating  statistically  identical  and  homogeneous  microstructures  with  smooth  gradient  in  pore 
size  resulting  from  changing  the  processing  parameters.  Strong  contrast  statistical  approach  is  also  used 
to  predict  the  in-plane  temperature  dependent  effective  electrical  conductivity  of  the  gradient  porous 
cathode  and  the  results  are  compared  to  the  experimental  data. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

SOFCs  are  highly  efficient,  environmentally  friendly,  and  fuel- 
flexible  electrochemical  devices  as  an  alternative  power  source 
for  next-generation  energy  systems.  Extensive  research  is  under¬ 
way  to  meet  the  goals  toward  commercialization  of  these  devices. 
One  major  goal  is  to  improve  the  performance  and  reliability  of 
the  cells  while  minimizing  costs  that  are  essentially  associated 
with  high  operating  temperatures  (~800  to  1000°C)  and  pro¬ 
cessing  of  SOFCs.  Previous  studies  have  shown  some  progress 
toward  this  goal  mainly  aimed  at  design  of  functionalized  materi¬ 
als/microstructures  that  operate  efficiently  at  lower  temperatures 
(<800  °C)  [  1  -4].  A  common  way  to  improve  the  functionality  of  indi¬ 
vidual  layers  at  low  temperature  is  to  apply  functionally  graded 
materials  (FGMs)  approach  with  compositional/microstructural 
variations  such  as  pore  size  in  some  spatial  direction  [5-9].  In 
composition  grading,  volumetric  ratios  of  the  electronic  conduc¬ 
tor  phase  is  mixed  with  ionic  conductor  electrolyte  to  form  the 
composite  electrode.  Furthermore,  particle/pore  grading  is  imple¬ 
mented  by  having  large  particles/pores  at  the  outer  layers  close  to 
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the  surface  for  gas  flow  and  gradually  smaller  pores  at  the  interface 
with  the  electrolyte  for  diffusion  and  chemical  reaction  processes. 
In  either  ways,  the  goal  is  to  enhance  the  electrochemical  reactiv¬ 
ity  and  increase  the  electronic/ionic  conductivity  of  the  electrode, 
to  minimize  the  mass  transport  resistance  to  gas  mixtures  and 
thus  concentration  polarization  inside  the  porous  electrode  and  to 
increase  mechanical  stabilities. 

Processing  of  the  SOFCs  electrodes  are  typically  carried  out  using 
conventional  processing  techniques  such  as  consolidation  and  sin¬ 
tering  of  ceramic  powder  [10,11].  However,  deposition  routes  take 
the  advantage  of  governing  the  processing  parameters;  thus,  one 
can  utilize  these  techniques  over  the  conventional  methods  to 
optimize  the  microstructure  of  the  SOFC  components  by  making 
desirable  variations  in  the  microstructure  [12,13].  Researchers  have 
attempted  to  identify  those  processes  that  have  a  significant  impact 
on  cell  performance  [1].  Here,  the  role  of  microstructure  is  very 
critical  as  it  acts  as  a  direct  link  between  processing,  properties  and 
performance  of  the  fuel  cell.  However,  all  current  manufacturing 
processes  are  stochastic  in  nature  [1  ];  thus,  an  effective  and  mean¬ 
ingful  understanding  of  the  processing-microstructure-property 
relationships  demands  establishment  of  a  framework  (consist¬ 
ing  of  development  of  manufacturing  processes  with  controllable 
processing  parameters)  that  can  account  for  these  variations  in 
constitutive  laws  [14]. 
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Fig.  1.  SEM  cross-section  image  of  the  functionally  gradient  porous  LSM  cathode 
deposited  on  top  ofYSZ  electrolyte.  Red  circles  show  presence  of  large  pores  on  the 
top  layer  and  gradually  smaller  pores  next  to  the  electrode/electrolyte  interface.  (For 
interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred 
to  the  web  version  of  this  paper.) 


To  this  end,  statistical  continuum  approaches  have  been  devel¬ 
oped  based  on  probability  functions  to  account  for  more  details  of 
microstructure  heterogeneities.  One  general  approach  for  predict¬ 
ing  the  effective  properties  of  a  two-phase  material  with  properties 
of  each  phase  similar  to  the  average  is  called  “weak-contrast” 
expansion.  However,  in  the  case  where  the  properties  of  the  phases 
are  highly  dissimilar,  “strong-contrast”  theory  is  implemented. 
While  most  of  the  analyses  reported  in  the  literature  are  performed 
on  finding  the  correlation  between  the  microstructure  and  prop¬ 
erties,  correlation  between  the  three  microstructure,  properties 
and  processing  parameters  is  less  regarded.  In  this  paper,  statis¬ 
tical  correlation  functions  are  used  to  make  a  linkage  between 
microstructure  of  a  functionally  graded  LSM  cathode  to  its  prop¬ 
erties  and  processing  conditions.  Gradient  porous  LSM  cathodes 
are  processed  by  a  modified  multi-step  spray  pyrolysis  [15].  Vari¬ 
ations  in  pore  geometry,  clustering  behavior  and  percolation  are 
also  studied.  Strong-contrast  statistical  continuum  theory  is  also 
used  to  predict  the  effective  conductivity  of  the  porous  LSM  cathode 
through  incorporation  of  the  details  of  microstructure  into  corre¬ 
lation  functions.  Results  are  then  compared  to  the  experimental 
measurement  data  of  the  electrical  conductivity. 


different  sizes  in  the  image;  note  that  the  SEM  image  is  taken  at  7° 
tilt  angle  for  clear  observation  of  the  gradient  in  the  pore  size  in  the 
cathode  layer. 


2.2.  Imaging  technique 

In  porous  brittle  ceramic  samples,  due  to  low  contrast  at  pore 
edges,  electron  microscopy  is  superior  to  optical  microscopy  meth¬ 
ods  that  involve  flat-polishing  step  for  quantitative  determination 
of  pore  dimensions  [16,17].  In  this  work,  the  in-plane  (, x-y  plane) 
microstructure  of  the  processed  gradient  porous  LSM  cathode  and 
its  temperature  dependent  electrical  conductivity  in  the  corre¬ 
sponding  direction  is  investigated.  For  this  purpose,  FEI  Nanolab 
Nova200  dual-beam  FIB/SEM  was  utilized  to  obtain  microscopic 
two-dimensional  (2D)  SEM  images  in  x-y  plane  by  sectioning 
the  gradient  porous  LSM  specimen  from  the  surface  down  to  the 
electrolyte  interface  in  vertical  direction  along  z  axis.  Using  Auto 
Slice  and  View  software  (FEI  Co.)  serial-sectioning  of  gradient 
porous  LSM  cathode  was  accomplished  to  observe  variations  in  the 
microstructure  in  x,  y,  z  directions.  The  dual-beam  FIB/SEM  is  com¬ 
posed  of  ion  beam  which  allows  milling  of  the  surface  while  the 
imaging  is  conducted  by  the  electron  gun.  A  number  of  LSM  serial- 
section  images  were  milled  at  highest  current  density  (20  nA)  and 
scanned  at  identical  spacing  of  100  nm  along  the  thickness,  i.e.  the 
gradient  (z)  direction.  Starting  from  the  rough  surface,  a  larger  area 
of  the  gradient  porous  LSM  sample  was  initially  milled  to  obtain 
a  flat  surface  for  further  accurate  milling  process.  A  set  of  seven 
2D-SEM  images  (50  pan  x  50  jxm)  at  equal  spacing  of  ~1.5  |xm  were 
chosen  among  serial-sections  for  characterization  and  as  represen¬ 
tative  of  the  variations  in  pore  size  in  the  microstructure.  These 
images  were  cut  out  of  a  larger  (>50  [xm  x  50  fxm)  area  to  prevent 
inappropriate  edge  effects  from  the  milling  process.  Then,  using  the 
threshold  technique,  images  were  converted  from  grayscale  into 
binary  (black-and-white);  black  pixels  are  associated  to  pore  phase 
with  a  value  of  zero  (phase  1 )  and  white  pixels  having  the  value  of 
one  represent  solid  LSM  phase  (phase  2).  Fig.  2  shows  the  schematic 
of  milling  process  of  gradient  porous  LSM  cathode  and  seven  binary 
images  obtained  for  further  microstructure  analysis.  As  shown,  the 
milling  process  is  carried  out  starting  from  surface  containing  large 
pores  and  ending  with  smaller  pores  at  the  electrolyte  interface. 

Every  two  binary  images  correspond  to  a  layer  that  has  been 
processed  at  a  specific  deposition  condition. 


3.  Microstructure  characterization  of  gradient  porous  LSM 
cathode 

3.1.  n-Point  correlation  functions 


2.  Experimental 

2.1.  Sample  preparation  procedure 

Continuously  graded  porous  cathode  functional  layer  was  pro¬ 
cessed  by  a  novel  multi-step  spray  pyrolysis  technique  reported 
elsewhere  [15].  The  LSM  was  deposited  on  a  1  mm  thick  pellet 
of  polycrystalline  YSZ  electrolyte  by  changing  the  deposition  vari¬ 
ables,  such  as  solution  flow  rate,  nozzle  to  substrate  distance  and 
temperature.  The  gradient  in  pore  size  is  achieved  by  design  of 
sequence  of  processing  steps  in  SP  including  variation  of  deposition 
parameters  in  a  controlled  manner  to  create  a  fine  microstructure 
containing  small  pores  and  nanocrystalline  particles  at  the  inter¬ 
face  and  gradually  larger  pores/particles  at  the  outer  layers  close 
to  the  surface  of  the  cathode  as  shown  in  Fig.  1.  Variation  in  pore 
size  in  the  gradient  microstructure  is  highlighted  by  red  circles  of 


Statistical  homogenization  techniques  use  n- point  correlation 
functions  to  statistically  describe  distribution,  orientation  and 
geometry  of  the  phases  in  the  microstructure  and  to  determine  the 
effective  physical  properties  of  the  heterogeneous  media.  In  the 
case  of  porous  media,  one  major  morphological  descriptor  of  the 
heterogeneous  microstructure  is  the  volume  fraction  of  the  phases 
(0)  or  porosity  that  can  be  determined  by  calculating  the  one-point 
correlation  function.  In  fact,  one-point  probability  function  repre¬ 
sents  only  the  volume  fraction  and  gives  no  more  information  about 
the  morphology  (Eq.  (1)).  Thus,  higher-order  probability  functions 
are  exploited  to  access  more  details  of  the  microstructure. 

P{(*)  =  0i  (1) 

Two-point  correlation  function,  P|(x,  x'),  is  another  probability 
function,  if  a  vector  r  is  thrown  on  a  two-phase  microstructure  such 
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Fig.  2.  (a)  Schematic  of  gradient  porous  LSM  cathode  on  YSZ  electrolyte,  (b)  Top  view  of  the  milled  area  of  the  gradient  porous  LSM  sample,  (c)  Set  of  binary  serial-section 
images. 


that  the  starting  point  of  this  vector  lies  at  point  x  in  phase  i  and  its 
end  lies  at  point  x!  in  phase  j: 

p|(x,  x!)  =  P  {(xg  i)  n  {x'  ej)}  i  =  1,2  and  j  =  1,2  (2) 

The  following  relations  hold  for  a  two-phase  microstructure: 

Pl2(x,  x ')  =  fa,  and  P%(x,  x')  =  (pj  ||r||  -»  0  (3) 

P%(x,x')  =  P\(x)xP{(x')  =  0i0j  ||r||  — ►  oo  (4) 

It  turns  out  from  (3)  that: 

P”+P;!2=01  (5) 

pf+pf  =  4,2  (6) 

P\2  =  Pf  (7) 

pp +P\2 +Pf +P?  =  1  (8) 

thus,  one  of  the  above  probabilities  comes  out  as  an  indepen¬ 
dent  variable.  These  correlation  functions  are  generally  analytically 
formulated  by  Corson  in  the  form  of  an  exponential  regression 
function  as  in  the  following  equation: 

P'i  (r)  =  4>i<t>j  +  (-P,+14>i<Pj  exp  (-c0rn« )  (9) 

where  Cy  and  riy  are  empirical  least  square  constants  of  the  curve 
fitted  to  the  experimental  probability  curves. 

3.2.  Pore-size  distribution  function 

The  main  concern  in  characterization  of  the  gradient  porous 
microstructure  is  to  realize  the  geometric  variations  of  the  pore 
phase  (such  as  size,  shape,  etc.)  along  the  gradient  direction  (z).  In 
general,  such  geometrical  information  is  not  contained  in  the  cor¬ 
relation  functions.  Based  on  a  definition  by  Torquato,  the  pore-size 
distribution  function  p(l)(r)  is  defined  by  the  probability  of  finding  a 
pore  with  its  length  between  r  and  r + dr  in  phase  (i)  of  interest.  This 
information  is  useful  for  studying  the  transport  properties  and  fluid 
flow  in  porous  media  [18].  Since  the  gradient  porous  microstructure 
is  statistically  isotropic  in  x  and  y  directions,  a  test  line  of  defined 
length  can  be  thrown  through  the  RVEs  in  either  isotropic  direc¬ 
tion  without  changing  the  orientation.  To  this  end,  a  MATLAB  code 
was  written  to  determine  the  length  of  the  phase  of  interest  (i.e., 
pores)  in  the  direction  of  the  test  line.  The  program  starts  with  scan¬ 
ning  the  RVEs  in  the  horizontal  (x)  direction.  The  number  of  line 
scans  is  equal  to  the  number  of  pixels  iny  direction  for  each  image. 


The  code  counts  the  number  of  pixels  in  the  given  line  segment 
for  its  binary  value  (denoted  by  0  or  1)  and  gives  the  number  of 
pixels  contained  completely  in  that  particular  phase.  By  normaliz¬ 
ing  the  values,  directional  pore-size  distribution  functions  for  given 
RVE  images  of  the  gradient  porous  sample  along  z  direction  can  be 
obtained. 

3.3.  Percolation  and  clustering  functions 

The  study  of  clustering  and  percolation  is  important  in  predic¬ 
tion  of  transport  properties  in  porous  structures.  It  is  performed 
by  introducing  a  two-point  correlation  cluster  function  C^\x  -  x') 
which  is  the  probability  that  both  points  x  and  x!  are  found  in 
the  same  cluster  of  one  of  the  phases  to  analyze  the  behavior  of 
conduction  and  flow  in  porous  microstructures.  The  information 
obtained  from  cluster  function  is  not  sufficiently  contained  in  any 
of  the  lower-order  microstructural  functions  described  above  [18]. 
Among  different  types  of  algorithms  used  for  solving  the  percola¬ 
tion  problem,  an  efficient  percolation  algorithm  is  used  to  check 
the  continuity  of  cells  for  all  clusters  [19]. 

3.4.  Strong  contrast  approximation  of  the  effective  conductivity 

Strong-contrast  approach  in  statistical  continuum  theory  is  used 
in  this  work  to  predict  the  effective  electrical  conductivity  of  the 
two-phase  isotropic  (porous)  media  based  on  the  work  of  Torquato 
[20],  where  the  individual  properties  of  phases  are  highly  dissim¬ 
ilar.  Mikdam  et  al.  have  reported  that  strong-contrast  formulation 
provides  a  better  estimate  for  the  effect  of  the  microstructure  on  the 
conductivity  for  porous  media  [21,22].  The  linear  constitutive  rela¬ 
tionship  can  be  generalized  to  describe  the  electrical  and  thermal 
conductivity. 

The  electrical  conductivity  a  at  any  arbitrary  point  x  is  related 
to  the  local  current  density  J{x)  and  the  local  electric  field  E(x)  by 
the  Ohm’s  law: 

i(x)  =  orc(x).£(x)  (10) 

The  local  polarization  field  P(x)  is  defined  by 
P(x)  =  [cre(x)-crR]E(x)  (11) 

Since  the  current  density  is  divergence  free,  we  will  have: 

V  •  J(x)  =  crRV  •  E(x)  +  V  •  P(x)  =  0 


(12) 
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Fig.  4.  One-point  correlation  function  (volume  fractions)  of  the  pore  phase  for  the 
microstructure  of  the  set  of  2D  SEM  micrographs  of  gradient  porous  LSM  cathode. 


4.  Property  measurements 


Fig.  3.  Representation  of  vectors  in  spherical  coordinate  [25]. 


4.1.  Electrical  conductivity 


Here,  a  potential  field  cp(x)  is  defined  such  that 

E(x)  =  V<p  (13) 

Thus,  by  substituting  Eq.  (13)  in  (12)  we  will  obtain  the  following 
equation: 

v  •  J(x)  =  -a*V  •  (</>(*))  +  V  .  P(x)  =  0  (14) 

The  solution  to  Eq.  (14)  can  be  obtained  using  Green’s  function 
and  integration  by  parts.  The  reader  is  referred  to  [23]  for  details 
on  derivation  of  the  solution.  Assuming  isotropic  properties  of 
the  phases,  the  effective  conductivity  tensor  (cre)  of  the  porous 
microstructure  in  any  d-dimensional  space  can  be  obtained  using 
the  strong-contrast  expansion  as  follows  [21,24]: 


Temperature  dependent  electrical  conductivity  of  the  half-cell 
LSM  cathode  was  measured  using  a  four-point  probe  method 
combined  with  a  micro-heater  for  heating  the  sample  up  to  its 
operating  temperature  (~700°C).  This  method  is  used  in  order  to 
obtain  the  in-plane  electrical  conductivity  of  the  processed  gradi¬ 
ent  porous  LSM  cathode  and  to  compare  the  predicted  in-plane 
values  for  corresponding  2D  SEM  images.  Measurement  was  per¬ 
formed  in  a  temperature  range  from  250  to  700  °C,  with  heating  and 
cooling  rate  of  10°Cmin_1  in  the  ambient  atmosphere.  The  auto¬ 
mechanical  stage  of  the  four-point  probe  setup  is  designed  such 
that  the  tips  can  move  up  and  down  to  make  sufficient  contact 
between  the  tips  and  the  porous  layer. 

5.  Results  and  discussion 


{<Te  —  Or/}  1  .{cre  +  (d  —  1  )(Tr/} 


1 


PsrP^x) 


I-aRd 


/ 


Pf(x,  x')  -  pS(x)pS(x') 
Pf(x)Pf(x  o 


HR{x  -  x')dx' 


-<JRd2pSR 


Pf(x,  x x") 
J*(x)i*(x') 


Pf(x,  x')P*{x',  x") 
pS(x)pS(x')P^(x") 


Hr(x  -  x') 


•  Hr(x'  -  x")  dx'  dx"  -  . . . 


(15) 


where  /  is  the  second-order  identity  tensor,  pSR  is  the  polarizability, 
crR  is  the  reference  conductivity,  and  HR(x-x')  is  a  second-order 
tensor,  defined  below.  The  polarizability  fiSR  is  expressed  by: 


Psr  = 


<?S  -  &R 

crs  +  (d-\)crR 


The  second-order  tensor  HR(x  -  x')  is  defined  by: 


(16) 


Hr(x  -  x')  =  — l- 

crR  \x  -  XT 


(17) 


where  is  the  total  solid  angle  contained  in  a  d-dimensional  sphere 
and 


n  = 


(*-*0 

\x-x'\ 


(18) 


Here,  the  three-point  correlation  function  is  estimated  using  the 
following  approximation  (with  the  representation  of  the  corre¬ 
sponding  vectors  in  spherical  coordinates  shown  in  Fig.  3)  [25]: 


PsJx,x',x")  = 


l**3 


\xx'\  +  \xx"\ 


PsJx,x")  + 


I**"! 


\XX'\  +  \XX"\  2 


PS2(X,X '] 


Ps2(x',x" 

P\(x) 


(19) 


The  plot  in  Fig.  4  shows  one  point  correlation  function,  i.e. 
variations  in  volume  fraction  of  the  pores  (porosity  %)  in  terms 
of  distance  from  the  surface  of  the  electrode.  Here,  the  average 
volume  fraction  of  the  pore  phase  is  0.54  which  is  obtained  by 
averaging  the  volume  fractions  of  the  over  the  individual  images. 
As  shown  in  this  plot,  the  volume  fraction  of  the  pores  is  not 
significantly  changing  through  the  microstructure  while  the  corre¬ 
sponding  SEM  micrograph  of  the  milled  area  reveals  that  the  pore 
size  is  decreasing  as  we  get  close  to  the  electrode/electrolyte  inter¬ 
face  in  z  direction.  Fig.  5  depicts  the  pore-size  distribution  (PSD)  for 
the  microstructure  of  the  set  of  2D  SEM  micrographs  of  the  gradi¬ 
ent  porous  LSM  cathode  shown  in  Fig.  2.  The  trend  shown  in  this 
plot  confirms  the  gradient  in  pore  size  in  the  microstructure  of  the 
fabricated  cathode  as  where  it  is  desirable  to  have  large  pores  to 
allow  mass  transport  from  the  outer  layers  to  inner  layers  with 
smaller  pores  at  the  electrode/electrolyte  interface  where  the  reac¬ 
tion  takes  place.  In  addition,  results  from  the  pore  size  distribution 
in  these  microstructures  reveal  that  the  number  of  pores  smaller 
than  1  |jim  is  increasing  from  surface  toward  the  interface  in  the 
direction  of  porosity  gradient,  thus  resulting  in  formation  of  more 
two-phase  boundaries  in  the  microstructure  for  electrochemical 
reactions.  Moreover,  it  is  clear  from  the  distribution  curves  that  the 
corresponding  phase  percolation  is  increasing  in  the  same  direction 
of  porosity  gradient. 

Fig.  6  shows  the  two-point  correlation  functions  (TPCF)  P22  or 
simply  P22  for  the  solid  phase  in  the  set  of  2D  SEM  micrographs 
shown  in  Fig.  2.  As  we  know,  when  r  approaches  to  zero,  the  prob¬ 
ability  function  gives  the  volume  fraction  of  the  corresponding 
phases.  The  TPCFs  shown  in  Fig.  6  are  overlapping  without  any  sig¬ 
nificant  variations  in  the  gradient  direction.  Results  show  that  all 
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Fig.  5.  Pore-size  distribution  functions  (PSD)  for  the  microstructure  of  the  set  of  2D 
SEM  micrographs  of  gradient  porous  LSM  cathode. 

images  have  statistically  similar  two-point  correlation  functions, 
morphologies,  and  volume  fractions,  except  for  the  size  distribu¬ 
tions.  Comparison  of  the  PSD  and  TPCF  results  also  indicates  that 
the  TPCFs  are  regardless  of  the  porosity  gradients  and  thus  the 
fabrication  method  (i.e.  SP)  is  capable  of  generating  homogeneous 
microstructures  that  are  statistically  identical  while  there  is  a  gra¬ 
dient  in  size  and  distribution  of  one  of  the  constituent  phases.  This 
also  confirms  anisotropy  in  z  direction  in  terms  of  pore  size  distri¬ 
bution  which  is  the  main  property  of  the  developed  gradient  porous 
cathode  microstructure. 

The  in-plane  clustering  of  the  solid  phase  C|2  is  shown  based  on 
two-point  correlation  cluster  function  (TPCCF)  in  Fig.  7.  As  shown 
in  this  plot,  at  higher  r  values,  the  most  isolated  in-plane  clusters 
attributes  to  Image  1  which  corresponds  to  the  top  surface  layer  of 
the  cathode  that  are  ineffective  in  conduction  while  percolation  in 
clusters  of  the  solid  phase  is  gradually  increasing  in  inner  layers  and 
is  observable  at  the  highest  values  in  Image  5.  It  should  be  noted 
that  the  results  presented  here  are  based  on  2D  analysis  of  the  in¬ 
plane  microstructure  of  the  cathode  and  due  to  columnar  structure 
of  the  processed  cathode,  through-plane  connectivity  of  both  solid 
and  pore  phases  are  assumed  to  be  preserved. 

Fig.  8a-c  shows  percolation  for  the  microstructure  of  the  set 
of  2D  SEM  micrographs  of  the  gradient  porous  LSM  cathode  in 
z  direction  from  surface  toward  the  electrode/electrolyte  inter¬ 
face.  This  figure  shows  clusters  of  solid  phases  in  different  colors 


Fig.  6.  Two-point  correlation  functions  (TPCF)  for  the  solid  phase  in  the  microstruc¬ 
ture  of  the  set  of  2D  SEM  micrographs  of  gradient  porous  LSM  cathode. 


Fig.  7.  Two-point  correlation  cluster  functions  (TPCCF)  for  the  solid  phase  in  the 
microstructure  of  the  set  of  2D  SEM  micrographs  of  gradient  porous  LSM  cathode. 

while  the  gradient  in  pore  size  as  well  as  the  particle  size  is  clearly 
observable.  The  percolated  phases  are  shown  in  the  same  color 
for  graphical  representation  and  observation  of  clustering  of  the 
solid  phase.  The  solid  phase  (LSM)  is  percolated  cluster  in  cho¬ 
sen  microstructures  in  the  porosity  gradient  direction  from  surface 
to  the  electrode/electrolyte  interface.  Connectedness  of  solid  pore 
phase  as  the  path  for  flowing  gas  is  shown  at  each  level  in  z 
direction  while  the  in-plane  pore  path  shrinks  and  distributes  into 
smaller  pores  at  inner  layers.  While  percolation  limit  is  about  30% 
which  defines  the  amount  of  porosity  required  for  a  3D  structure  to 
achieve  cathode  requirements  and  properties  at  the  optimum  level 
and  also  for  the  gas  diffusion  to  be  not  rate  limiting,  this  value  is 
higher  for  a  2D  structure  [26,27]. 

While  high  electrical  conductivity  of  the  LSM  cathode  at 
high  temperatures  results  in  lower  ohmic  polarization,  sufficient 
amount  of  porosity  in  the  microstructure  is  required  to  provide  a 
path  for  mass  transfer  through  the  cathode  microstructure.  Thus, 
it  is  important  to  know  how  the  morphological  changes  in  the 
microstructure  affect  the  electrical  conductivity.  In  order  to  predict 
the  temperature  dependent  electrical  conductivity  of  the  gradient 
porous  LSM  cathode,  electrical  conductivity  of  the  actual  LSM  phase 
with  the  same  chemical  composition  of  Lao.8Sr0.2Mn03+5  as  the 
processed  LSM  is  taken  at  different  temperatures  from  the  litera¬ 
ture  as  the  input  data  for  the  model  [28].  Fig.  9  shows  comparison  of 
the  predicted  electrical  conductivities  and  the  experimentally  mea¬ 
sured  values  as  a  function  of  reciprocal  temperature  in  air  for  the 
gradient  porous  LSM  cathode  after  heat  treatment  in  the  tempera¬ 
ture  range  of 200-700  °C.  As  shown,  the  electrical  conductivities  are 
increasing  by  increasing  the  temperature.  The  expression  for  the 
temperature  dependent  electrical  conductivity  of  the  small  polaron 
materials  is: 

where  A  is  the  pre-exponential  factor,  T  is  the  absolute  tem¬ 
perature,  I<  is  the  Boltzmann  constant,  and  Ea  is  the  hopping 
activation  energy.  In  general,  the  activation  energy  for  electron 
hopping  depends  on  the  temperature  and  material  structure.  In  the 
La0.8Sr0.2lVInO3  system,  the  substitution  of  Sr2+  in  La3+  results  in  the 
formation  of  Mn4+  ions,  some  of  which  in  combination  with  Mn3+ 
ions  form  Mn3+-0-Mn4+  bonds  that  mostly  contribute  in  electronic 
conduction.  As  seen  in  Fig.  9,  the  slopes  of  the  two  lines  repre¬ 
senting  the  activation  energy  for  electrical  conduction  match  in 
experimental  and  modeled  data.  However,  the  slope  of  the  mea¬ 
sured  data  is  slightly  increasing  at  400  °C  which  can  be  attributed 
to  the  change  of  conduction  mechanism  due  to  change  in  chemi¬ 
cal  state  of  the  material  such  as  formation  of  oxygen  vacancies  and 
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Fig.  8.  Percolation  for  the  microstructure  of  the  set  of  2D  SEM  micrographs  of  the  gradient  porous  LSM  cathode  changing  in  the  thickness  direction  from  surface  to  the 
electrode/electrolyte  interface  in  (a)-(c). 


their  contribution  in  ion  conduction.  The  close  values  of  the  acti¬ 
vation  energies  verify  similarity  in  stoichiometric  composition  of 
the  processed  LSM  material  and  confirm  that  the  proposed  model 
can  only  account  for  the  effect  of  microstructural  changes  in  pre¬ 
dicted  conductivity  behavior  [29].  Moreover,  observation  of  a  linear 
behavior  in  both  plots  reveals  that  the  hopping  mechanism  of  small 
polarons  between  M3+  -  M4+  pairs  in  this  temperature  range  is  the 
conduction  mechanism  in  p-type  perovskite  oxide  structure  LSM. 
However,  from  the  microstructural  point  of  view,  it  is  observed  that 
the  temperature  dependent  conductivities  of  the  gradient  porous 
LSM  cathode  are  slightly  lower  than  the  predicted  values;  con¬ 
sidering  the  amount  of  porosity  in  the  processed  cathode,  these 
values  are  still  in  reasonable  agreement  with  various  experimen¬ 
tal  data  reported  in  the  literature  and  required  for  SOFC  cathode 
application  [30,31].  Due  to  unavailability  of  data  on  temperature 
dependent  electrical  conductivity  for  the  exact  composition  of  the 
processed  sample,  the  input  values  taken  from  the  literature  may 
not  exactly  represent  the  actual  properties  of  the  solid  phase  with 
different  morphology.  Moreover,  percolation  is  assumed  to  be  a 
source  of  discrepancy  between  the  experimental  and  theoretical 
data  in  direction  of  conductivity.  Generally,  in  single-phase  perco¬ 
lation,  since  the  percolation  limit  for  a  2D  square  lattice  is  higher 
than  that  in  a  3D  cubic  one,  the  percolation  in  3D  is  more  read¬ 
ily  achievable  than  in  two  dimensions.  Here,  due  to  anisotropy 
in  the  microstructure  in  z  direction,  the  actual  percolation  limit 
which  is  strongly  dependent  on  the  distribution  of  the  conductive 
phase  in  specific  direction  is  different.  However,  due  to  branch-like 
(columnar)  morphology  of  the  processed  gradient  porous  cathode 
and  formation  of  some  gaps  between  these  columns,  the  in-plane 
2D  analysis  of  the  gradient  porous  microstructure  investigated  in 
this  work  as  well  as  different  percolation  limit  of  the  conductive 
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Fig.  9.  Comparison  between  the  experimental  and  predicted  values  of  the  temper¬ 
ature  dependent  electrical  conductivity  of  gradient  porous  LSM  cathode. 


phase  in  the  in-plane  direction  rather  than  conductance  direction 
have  resulted  in  overestimation  of  the  electrical  conductivity  by  up 
to  an  order  of  magnitude.  Therefore,  the  model  seems  to  predict 
a  different  conductivity  for  the  gradient  porous  LSM  than  what  is 
actually  measured.  In  a  recent  work  by  Martinez  et  al.,  the  perco¬ 
lation  physics  and  a  similar  effect  on  conductivity  behavior  is  also 
investigated  [32]. 

6.  Conclusion 

The  statistical  continuum  mechanics  is  used  to  correlate  prop¬ 
erties  of  the  gradient  porous  LSM  cathode  to  its  microstructure. 
Statistical  characterization  of  the  microstructure  revealed  the  pres¬ 
ence  of  pore  size  gradient  in  the  microstructure,  starting  from  large 
pores  at  the  electrode  surface  and  gradually  decreasing  the  pore 
size  toward  the  electrode/electrolyte  interface.  Results  of  statis¬ 
tical  characterization  also  show  that  all  images  have  statistically 
similar  two-point  correlation  functions,  morphologies,  and  volume 
fractions,  except  for  the  size  distributions.  Comparison  of  the  PSD 
and  TPCF  results  also  indicates  that  the  TPCFs  are  regardless  of  the 
porosity  gradients  and  thus  the  applied  fabrication  method  (i.e. 
multi-step  spray  pyrolysis)  is  capable  of  generating  homogeneous 
microstructures  that  are  statistically  identical  while  they  contain  a 
gradient  in  size  and  distribution  of  one  of  the  constituent  phases. 
Results  from  TPCCF  also  reveal  gradient  in  cluster  size  solid  phase 
as  well  as  the  pore  phase  in  z  direction.  Experimentally  measured 
temperature  dependent  electrical  conductivity  data  show  slightly 
lower  values  than  the  predicted  ones.  This  discrepancy  can  be 
referred  to  the  effect  of  morphology  and  percolation  effects  in  the 
porous  microstructure  as  well  as  unavailability  of  exact  conduc¬ 
tivity  data  for  the  composition  of  interest.  Columnar  structure  of 
highly  porous  gradient  microstructure  is  also  suggested  as  another 
fact  responsible  for  the  lower  values  of  the  conductivity  in  in-plane 
direction. 
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